Fundamental Quantum Limit to Waveform Estimation 
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We derive a quantum Cramer-Rao bound (QCRB) on the error of estimating a time-changing 
signal. The QCRB provides a fundamental limit to the performance of general quantum sensors, 
such as gravitational-wave detectors, force sensors, and atomic magnetometers. We apply the QCRB 
to the problem of force estimation via continuous monitoring of the position of a harmonic oscillator, 
in which case the QCRB takes the form of a spectral uncertainty principle. The bound on the force- 
estimation error can be achieved by implementing quantum noise cancellation in the experimental 
setup and applying smoothing to the observations. 
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The accuracy of any sensor is limited by noise. To 
quantify the potential performance of a sensor, it is often 
useful to compute a lower bound to the error in the esti- 
mation of the signal of interest. One of the most widely 
used bounds is the Cramer-Rao bound (CRB), which lim- 
its the mean-square error in parameter estimation [l[ . 

The development of quantum technology highlights the 
question of how quantum mechanics impacts the per- 
formance of sensors. Helstrom formulated a quantum 
Cramer-Rao bound (QCRB) which stipulates that 
the minimum estimation error is inversely proportional 
to a property of the sensor known as the quantum Fisher 
information. The QCRB is central to quantum sen- 
sor design in the burgeoning field of quantum metrol- 
ogy S 3 for several reasons. It allows one to determine 
whether the fundamental sensitivity of a sensor design 
meets the requirements of an application, provides a cri- 
terion against which the optimality of quantum sensing 
schemes can be tested, and motivates improvements of 
schemes that are suboptimal. For sensors near the fun- 
damental limit, the QCRB can also be used to quantify 
the trade-off between sensing accuracy and physical re- 
sources of the sensor, so that efficient ways of improving 
sensitivity can be identified. 

Most prior work on the QCRB considered estimation 
of one or a few fixed parameters. Yet, in most sensing ap- 
plications, such as force sensing and magnetometry, the 
signal of interest is changing in time. This time-changing 
signal, which we call a waveform, is coupled continuously 
to the sensor, and continuous measurements on the sensor 
are used to extract information about the waveform 0- 
0]. Here we derive the QCRB for waveform estimation — 
the first such derivation to our knowledge — allowing for 
any quantum measurement protocol, including sequen- 
tial, discrete or continuous measurements. 

Previous work on the QCRB generally did not take 
into account prior information, but for the task of esti- 
mating a waveform, which often depends on an infinite 



number of unknown parameters, parameter estimation 
techniques no longer suffice and prior information is re- 
quired to make the problem well defined The prior 
information might, for example, restrict the signal to a 
finite bandwidth, making integrals over frequency finite 
that otherwise would diverge. Thus a crucial feature of 
our QCRB is the inclusion of prior waveform information. 

Our result provides a rigorous criterion against which 
the optimality of design, control, and estimation strate- 
gies for quantum sensors, such as gravitational-wave de- 
tectors, force sensors, and atomic magnetometers, can be 
tested. As an example, we calculate the QCRB on the 
error of force estimation via continuous position measure- 
ments of a harmonic oscillator, in which case the bound 
takes the form of a spectral uncertainty principle. We 
show that the bound can be achieved by implementing 
quantum noise cancellation (QNC) to remove the back- 
action noise from the observations Q and applying the 
estimation technique of quantum smoothing [7[ to the ob- 
servations. This proves the optimality of such control and 
estimation techniques for force sensing and establishes 
our QCRB as the fundamental limit to force sensing. 

Let x{t) denote the classical waveform to be estimated. 
For simplicity, we assume x{t) to be a scalar function; 
generalization to multiple processes is straightforward. 
We discrctize time as tj = to + J<5t, j = 0, 1, . . . , J, and 
assume that St is small enough that we can treat x{t) 
as picccwise-constant, i.e., x(t) = Xj for tj < t < tj+i- 
The prior probability density P[x] for the vector x = 
{xj-i, . . . , xq)'^ characterizes what is known or assumed 
about the waveform prior to the measurements. For a 
vector of observations y = {y^-i, ■ ■ ■ , yi,yo)'^ made any 
time during the interval Iq < t < tj, we define a condi- 
tional probability density The joint probability 
density is P[y,x] = P[?/|x]P[a:;]. Finally, we define the es- 
timate of Xj as Xj [y] and the estimate bias, given signal x, 
as jDy {xj - Xj)P[y\x] = bj[x], where Dy = JJ^'q dyn- 

Multiplying both sides of hj[x\ by P[x\, differentiating 
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with respect to a;^, and then integrating over ah x using 



Dx = n^=o^ '^^j-i obtain 



J Dx Dy {xj — Xj 



dP[y, 



dxk 



= lDx-^{b,[x]P[x])=0, (1) 

where the final cquahty assumes bj[x]P[x]\x^=±oo = 0. 
This assumption, also used in the proof of the classi- 
cal CRB is satisfied as long as the prior density ap- 
proaches zero at the infinite endpoints (as it must for any 
probability density) and the bias there is not infinite. 

Quantum mechanics enters this description, which till 
now is classical, by determining the conditional probabil- 
ity of the observations. Given a quantum system, we can 
describe any measurement protocol, including sequential 
measurements and excess decoherence, during the inter- 
val to < t < tj hy introducing appropriate ancillae, in 
accord with the Kraus representation theorem 
This also accounts for any feedback during the interval, 
based on the measurement outcomes, because the princi- 
ple of deferred measurement allows one to put off the 
measurements on the ancillae till time tj] measurement- 
based feedback is replaced by controlled unitaries prior 
to the measurements, as schematically shown in Fig. [TJ 
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FIG. 1: Any quantum dynamics and sequential measurements 
described by completely positive (CP) maps, including feed- 
back based on measurement outcomes, as illustrated in (a), 
can be reproduced by unitary evolution of an enlarged system 
that includes appropriate ancillae, coherent controlled uni- 
taries, and deferred measurements of the ancillae, as shown 
in (b). 

In this approach, the overall system dynamics is de- 
scribed by unitary evolution of the enlarged system; 
the conditional probability of observations is given by 
P[y|a;] = tr(£'[y]p^), where px is the density operator of 
the enlarged system at time ij, conditioned upon x, and 
E[y] is the positive-operator-valued measure (POVM) 
that describes the (deferred) measurements up to time tj. 
We denote expectation values with respect to px by angle 
brackets subscripted by x, so that {E[y])x = tr(^E[y]px) ■ 
Continuous measurements can be modeled as the limit of 
a sequence of infinitesimally weak measurements Q. 

We now follow a procedure similar to the one used 
by Helstrom to derive the QCRB. We introduce an 
operator Qk that satisfies dpx/dxk = {QkPx + PxQ\)/'^- 
Unlike Helstrom, we do not require Qfc to be Hermitian. 



Note that the vanishing trace of dpx/dxk in the definition 
of Qk implies that Rc{Qk)x = 0. 

It is convenient to incorporate the prior information 
by working in terms of a density operator p[x] = pxP[x] 
in a hybrid quantum-classical space and introducing 
an operator Lk = Qk + dlnP[x]/dxk, which satisfies 
dp[x]/dxk = {Lk[x]p[x] + p[x]lI[x])/2. In terms of Lfe, 
Eq. Il]) takes the form that we use to derive the QCRB: 



Sjk = Re J DxDy {xj - Xj) tT(E[y]Lk[x]p[x]) 



(2) 



Multiplying Eq. ^ by UjVk, where Uj and Vk are the 
components of arbitrary real column vectors u and v, 
and then summing over all j and fc, we obtain 



T 
V u 



^UjVj = Re j DxDy tv{A^ B) , 



(3) 



where A"^ = Y.k y^E[y]Lky/p\x], B = Y^^Ujixj - 
Xj')^ p\x^^ E\y\ and T denotes transposition. It follows 
from Eq. ^ that 



(v^uf < J DxDy tr(AtS) 

< J Dx Dy tiiA'' A) J Dx Dy tiiB^ B) , (4) 

where the second inequality is the Schwarz inequality. 

The second integral in Eq. (j4|) is j Dx Dy tr{B^ B) = 
u^'S u, where 



Ejfc = J Dx Dy P[x,y]{xj - Xj){xk - Xk) (5) 

is the estimation-error covariance matrix. The first inte- 
gral in Eq. (j4]) is, using the completeness of the POVM, 
jDxDy ti{A''A) = v'^Fv, where f is a (real, symmet- 
ric) Fisher-information matrix. 



F,k = \j DxP[x\ iv{{L]Lk + lILj)p, 



(6) 



Since Rc{Qk)x = 0, separates neatly into a quantum 
and a classical, prior-information component, i.e., F = 
F(Q) +F^^\ where 



F, 



AQ) 
jk 

(C) 



jk 



DxP[x] triiQlQk 



DxP[x] 



QIQ,)P.] (7) 
d\nP[x] d\nP[x] 



dxj dxk 



(8) 



When these results are substituted into Eq. (|4]) , we find 
that {v^ Fv){u^Y.u) > {v'^u){u^v). Setting v = P-'^u 
implies that v?" {Yj — F^^)u > for arbitrary real vectors 
u. Since S — F~^ is real and symmetric, this implies that 
I] — F~^ is positivc-semidefinite; the matrix inequality 



j:> p- 



(9) 
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is the QCRB in its most general form. To use a CRB 
in practice, it is customary to define a non-negative, 
quadratic cost function C = tr(A-^S) using a positive- 
semidefinite (Hermitian) cost matrix A suited to the ap- 
plication [l|, Q. The matrix QCRB is equivalent to a 
lower bound, C > tr(A'^i^^^), on aU such cost functions. 

To calculate the QCRB, we must be more specific 
about the evolution of the enlarged quantum system. 
The Hamiltonian that governs overall system dynam- 
ics over the interval tj < t < tj+i, of duration St, 
is Hj(xj), with corresponding evolution operator Uj ~ 
cxp[—iHj{xj)5t/h]. We have dUj/dxj ~ Uj{—ihj 6t/h), 
where hj = dHj/dxj. Let Ukj = Uk-i ■ ■ ■ Uj denote the 
evolution operator over the interval tj < t < tk- The 
density operator px is related to the initial density op- 
erator po by Pa; = UjoPqUIq, which gives dp^/dxh = 
-i[Mk,Px], where 



dxk 



St. 



(10) 



with hk 



uLhkUko 



Ik = u^^nkUkQ = h(tk) being the Heisenberg- 
picturc version of hk- An obvious choice for is the 
anti- Hermitian Qk — — 2iAA/fc, where AMk = Mk — 
{Mk)x- The quantum part of the Fisher matrix then be- 
comes 



jk 



^^"^^^ fDxP[x]^tT(^{AhjAhk + AhkAh^)po 



(11) 



where Ahk = hk — (hk)o- Angle brackets with subscript 
denote an expectation value with respect to po. The 
quantum Fisher information is thus a two-time covari- 
ance function, averaged over P[x]. 

To take the continuous-time limit, we let St 0, 
S.fc ^ F,k/{Stf F{t,,tk), and Ajk/iSt)^ ^ 

A{tj,tk)- The estimation-error covariance matrix be- 
comes the two-time covariance function of estimation er- 
ror, T,{t,t'), and the Fisher matrix becomes F{t,t') = 
F^Q\t,t') + F^^^{t,t'), with 

F^Q\t,t') = ^ JdxP[x] 

x-(^Ah{t)Ah{t') + Ah{t')Ah{t)^ , (12) 

^/--W^i^^i^. ° (13, 

5 / 5x{t) being the functional derivative. 

In the continuous-time limit, the matrix QCRB retains 
the same form as Eq. ([5]), where the continuous-time in- 
verse is defined by j/^' dt" F{t,t")F-\t" ,t') = S{t - t'). 
The bound on a cost function becomes 



C= / dtdt' A{t,t')J:{t,t') > I dtdt' A{t,t')F~\t,t') . 

(14) 



Equation ([M]), valid for any cost function, is the most 
serviceable expression of our chief result. An important 
special case is the point estimation error, 

n(t) = s(t,t) = ([^(^) - x{t)]^) > F~^{t,t) , (15) 

where angle brackets without a subscript denote an over- 
all quantum-classical average. 

To illustrate the use of our QCRB, we consider the esti- 
mation of a force x{t) on a quantum harmonic oscillator. 
The Hamiltonian is H = /2m -\- muj'j^q^ /2 — qx{t), with 
q being the position operator, p the momentum operator, 
m the mass, and ujm the resonant frequency. In this situ- 
ation, we have h{t) = dH{x{t)) / dx{t) = — g, which leads 
to a quantum component of the Fisher information. 



F^^\t,t') = ~(^Aq{t)Aq{t') + Aq{t')Aq{t) 



(16) 



The further average over P[x] in Eq. (|11|) can be omitted 
in Eq. (|16p because x{t) appears linearly in q{t) and thus 
drops out of Aq{t). If we assume that x{t) is a Gaussian 
process, the classical, prior-information component of the 
Fisher information is the inverse of the prior two-time 
covariance function of Ax(t) [l|. 

We now assume that all noise processes are stationary. 
For a stationary, zero- mean process /(t), the covariance 
function {f{t)f{t')) depends only on the time difference 
T = t' ~ t and can be Fourier-transformed to give the 
power spectral density S'/(a;) = J^^dT{f{t)f{t+T))e^'^^. 
The choice A{t,t') = exp[ia;(i' - t)]/{tj - to), together 
with taking to ^ —oo and tj — > oo, makes C{uj) the 
power spectral density of the estimation error. The 
QCRB (|14p then becomes a spectral uncertainty princi- 
ple: 



C(w)h5A,(^) + 



> 



(17) 



In the time-stationary case, the matrix QCRB is equiv- 
alent to satisfying this spectral uncertainty principle for 
all cj. A bound on the point estimation error now follows 
fromH = /^(dw/27r)C(w). 

To proceed in our approach, we must specify the mea- 
surements that extract the force information from the 
oscillator and include the associated backaction. Thus 
we now suppose that one performs continuous position 
measurements, using, for example, a continuous opti- 
cal probe. The observation process is y = q + rj, and 
the oscillator equations of motion are dq/dt = p/m and 
dp/dt = —muj'^q + x + S,, where ^ is the backaction noise. 
Here £^{t) and r](t) are like the quadrature components of 
an optical field, obeying the canonical commutation re- 
lation [^{t), ri{t')] = ihS{t — t'). We assume ^ and 77 have 
zero mean; their spectra satisfy an uncertainty principle. 

If we introduce a small amount of damping, Aq{t) 
is the inhomogeneous solution for q(t), driven just by 
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£,{t), and becomes stationary. In the limit of negligi- 
ble damping, the spectrum of Aq becomes S'Ag(w) = 
\G{u})\'^S^(u!), where G(w) = l/m(a;,^„ - w^) is the os- 
cillator transfer function. The spectral uncertainty prin- 
ciple pT)) now takes the form 



> 



45Ax(c^)y - 4 

The corresponding bound on point estimation error is 



(18) 



n > 



^oo 27r V« OAx(w) 



(19) 



Notice that a bandwidth constraint on x(t) is incorpo- 
rated in the prior information: Sax{^) goes to zero out- 
side the relevant bandwidth, thus allowing C{uj) to be 
zero there and making the integral ([T9l) finite. 

We can elucidate the meaning of the QCRB ([T^ by 
considering how to estimate the force from the obser- 
vations in this scenario. In the frequency domain, the 
observation process y{t) reads y{u}) = G{uj)[x{uj) + z{uj)], 
z being a noise term that depends on ^ and r/. Using 
smoothing P, Q to estimate x from y yields an error 



n= 

Loo 2^ 



(20) 



This is the minimum achievable error for a given noise 
spectrum Sz(j-o)- It cannot be reached by the more well- 
known technique of filtering Q , as filtering does not make 
use of the entire observation record. If and rj arc un- 
correlated and quantum limited, we have 



|G(c.)P 



\GH\ 



= Ssql{co), (21) 



where the power spectrum S'sql('^) is known as the stan- 
dard quantum limit (SQL) for force detection 

It is now evident that to attain the QCRB (fT9|) . it 
is necessary to beat the SQL. This requires evading 
or tempering the effects of the backaction ^. One way 
to do this is to correlate ^ and 77, as was proposed 
for interferometric gravitational-wave detectors by Un- 
ruh [ll| . An alternative is to use quantum noise cancel- 
lation (QNC) 0, which has the advantage of making the 
QCRB pop achievable, as we now show. One QNC ap- 
proach, discussed in adds an auxiliary oscillator with 
position q' and momentum p' . One monitors continu- 
ously the collective position Q ^ q + q' , giving a process 
observable y = Q + rj; the backaction force ^ acts on 
P = {p + p')/2 and thus equally, with strength ^, on each 
of the two oscillators. Suppose the auxiliary oscillator 
has the same resonant frequency and equal, but opposite 
mass (the negative mass can be simulated by an opti- 
cal mode at the red sideband of the optical probe). The 
dynamics of the collective position is then determined 
by dQ/dt — 6p/m and dSp/dt = —mio^^Q + x, where 



5p = p — p' . There being no backaction noise in z{t), one 
easily finds that 



> 



1 



>5',)(^) 

|G(c.)P - 45^^) \G{u)\ 



(22) 



with equality for quantum-limited noise. This quantum- 
noise-cancellation scheme beats the SQL and if the noise 
is quantum limited, docs so optimally: the smoothing 
error given by Eq. (PU)) achieves the QCRB p^ . which 
implies that the spectral uncertainty principle (jlSp is sat- 
urated. Our force-sensing QCRB, rigorously proven and 
demonstrably achievable, thus serves as a fundamental 
quantum limit, against which the optimality of future 
force sensing schemes should be tested. More generally, 
our QCRB for arbitrary cost functions will find ap- 
plication whenever quantum-limited estimation of tem- 
porally varying waveforms is attempted. 
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